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The electronic structure of (group-III) metal-monochalcogenide monolayers exhibits many un¬ 
usual features. Some, such as the unusually distorted upper valence band dispersion we describe as 
a ‘caldera’, are primarily the result of purely orbital interactions. Others, including spin splitting 
and wavefunction spin-mixing, are directly driven by spin-orbit coupling. We employ elementary 
group theory to explain the origins of these properties, and use a tight-binding model to calculate the 
phenomena enabled by them, such as band-edge carrier effective g-factors, optical absorption spec¬ 
trum, conduction electron spin orientation, and a relaxation-induced upper valence band population 
inversion and spin polarization mechanism. 


I. INTRODUCTION 

Three-dimensional van der Waals solids from the 
(group-III) metal monochalcogenide layered semiconduc¬ 
tor family MX (where M is Ga or In, and X is S, Se, 
or Te) have been intensively investigated by both exper¬ 
iment and theory for many decades. The classical litera¬ 
ture on this subjectP contains reports of measurements on 
optical absorption, photo- and electroluminescence, pho¬ 
toconductivity, radiative recombination, electrical con¬ 
ductivity, and Hall effedP^^. Even conduction electron 
spin polarization via optical orientation w as car ried out 
to study spin-dependent carrier dynamics. ^^dd ^he 
theory side, we find a first attempt to derive the band 
structure based on symmetry analysis and a simple few- 
band model occurring over fifty years ago.^ Since then, 
many detailed bandstructure calculations have appeared, 
utilizing the empirical pseudopotential methocP^®^ or 

tight binding formalism.^SHM] 

However, it was not until the recent search for 
beyond-graphen^^Hm two-dimensional semiconductors 
(such as transition metal dichalcogenide^^lHlIl and 
phosphorend22H2^ that this class of material was exper¬ 
imentally realized down to few- or mono-layer thickness 
by mechanical exfoliation from bulk crystals. Experimen¬ 
tally, GaS and GaSe ultrathin layer transistors have been 
demonstrated.!^ Photoluminescence measurements show 
exciton features and reduc tion o f optical efficiency when 
sample thickness decreasesGircularly polarized pho¬ 
toluminescence reveals spin dynamics in nanoslabs. 

In some cases, monolayer MX can even be synthesized 
epitaxially on silicoii^ or non-epitaxially on insulating 
substrates such as Si 02 via vapor pha se de position, with 
quality that rivals exfoliated material .1^2111] 

Despite this recent explosion of experimental results 
with monolayer MX, the theoretical establishment has 
relied almost exclusively on sophisticated ab initio meth¬ 
ods to model electronic structure. In particular, an un¬ 
usual distortion of the highest valence band (sometimes 
called an ‘upside-down Mexican hat^^ or ‘sombrero’) is 
predicted to create an indirect bandgap and a density-of- 
states singularity at or near the band edge.!^ This fea¬ 
ture, which is more appropriately called a ‘caldera’, van¬ 


ishes in the bulk.^^l First-principles schemes have made 
other fascinating predictions, suc h as s pontaneous mag¬ 
netism in p-type monolayer GaSe.!^^!^ However, the un¬ 
derlying fundamental physics at the root of many intrigu¬ 
ing properties is obscured by these brute-force numerical 
approaches, leaving many elementary questions without 
satisfactory answers. 

Our aim in the present work is to investigate the under¬ 
lying symmetries responsible for the many unique prop¬ 
erties that are common to all MX monochalcogenide 
monolayers sharing the same type of lattice structure and 
zone-center band-edge states. Using elementary group 
theory, we reveal the origin of the extraordinary ‘caldera’ 
shape valence band edge, and examine various important 
phenomena made possible by the spin-orbit symmetry. 
These include orbital degeneracy breaking in the valence 
band, fc-cubic dependence of the lowest order Dressel- 
haus splitting, orbital magnetism and effective Lande g- 
factor, and spin dynamics during optical orientation. To 
assist the reader in acquiring an intuitive and quanti¬ 
tative understanding of our theory, we take the case of 
GaSe as a specific example in numerical calculations, us¬ 
ing an empirical tight-binding model following Refs. [22] 
and |23] and incorporating on-site spin-orbit coupling 
parameters,!^ as well as ab initio density-functional the¬ 
ory (DFT) with the Quantum ESPRESSO package.!^ 
We emphasize that these numerical procedures are im¬ 
plemented only as a verification of the underlying physics 
determined by symmetry, which remains robust regard¬ 
less of numerical details (such as the choice of functional 
or pseudopotential in DFT). 

This paper is organized as follows: We start with back¬ 
ground information on the essential group theory in Sec. 
im and with the assistance of the nearly free electron 
model, we analyze the symmetry of the spin-independent 
bandstructure, revealing the interactions causing an un¬ 
usual valence band distortion: the ‘caldera’. In Sec. m 
we include spin-orbit interaction and investigate its effect 
on broken band degeneracy, eigenfunction composition, 
and spin splitting. In Sec. |IV| we show how orbital dia¬ 
magnetism and valley-spin coupling result from the spin- 
mixed wavefunction symmetry. Finally, in Sec. jV] we 
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discuss conduction electron spin orientation by optical 
excitation, including relaxation dynamics causing spin- 
polarized population inversion in the lower valence band 
and three-level spin pumping of the upper valence band. 

II. SYMMETRY SANS SPIN 

The unit cell of the monolayer metal- 
mono chalcogenide hexagonal lattice is composed of 
two group-III metal and two chalcogen atoms, forming 
an upper and lower sublayer related by in-plane mirror 
reflection symmetry [Fig. 0a)]. Within each layer, the 
two types of atoms are covalently bonded and arranged 
alternately at the honeycomb lattice sites [Fig. 0b)]. 
The honeycomb lattice is buckled such that metal atoms 
are closer to the opposite sublayer. In this configuration, 
the two sublayers are tightly bound by adjacent metal 
atoms [blue bond in Fig. 

The first Brillouin zone is shown by the hexagon in¬ 
scribed within the reciprocal lattice in Fig. 0c), together 
with high symmetry points F, K {K') and M. Reciprocal 
lattice points equidistant from the origin can be divided 
into three sets according to their symmetry, denoted in 
Fig. 0 c) by three different types of markers (green hex¬ 
agram, blue square and red circle), corresponding to the 
three types of zone-center states in the nearly free elec¬ 
tron band structure shown in Fig. 0d). Before proceed¬ 
ing with a detailed discussion of electronic properties, 
we first present a brief symmetry analysis of the system 
using group theory. 

A. Group theory 

The space group of single-layer MX is symmorphic, 
with its point group including twelve symmetry op¬ 
erators divided into six classes, denoted by Ci...,6 in Ta¬ 
ble |T] The identity operator E is in class Ci. C2 includes 
the ±120° rotations along an out-of-plane axis [through 
the position of any atom or the honeycomb center in Fig. 
0b)]. C3 is composed of a 180 ° rotation along the three 
axes [solid lines in Fig. 0b)] within the plane bisecting 
the two sublayers. C4 includes the in-plane mirror reflec¬ 
tion (T/j. C5 and Ce take into account operators from the 
product of (Jh and those in C2 and C3, respectively. 

The Brillouin zone center F-point has the same sym¬ 
metry of the point group D^h- The six irreducible repre¬ 
sentations (IRs), denoted by r]*"2 3) have the characters 
given in Table |T] The plus (minus) superscript reflects 
the even (odd) parity with respect to uu- Correspond¬ 
ing to each IR, we give in Table the lowest order basis 
functions which we will use to describe the F-point states 
based on their wavefunction symmetries. 

For intuitive understanding, and to aid numerical cal¬ 
culation of the band structure in the tight-binding for¬ 
malism, we also provide in Table |l] the symmetries of s 
and p atomic orbital combinations, following Ref. [^ . 



FIG. 1 . Unit cell of single-layer group-III metal monochalco- 
genide in (a) perspective- and (b) plan-view. The red and 
green spheres correspond to chalcogen anions and metal 
cations, respectively. The dashed blue frame represents the 
unit cell boundaries. The three solid lines in (b) are axes 
for 180 ° in-plane rotation operations. Panel (c) shows re¬ 
ciprocal lattice points and reduced Brillouin zone, where dif¬ 
ferent colored markers indicate symmetry-related points of 
different zone-center character, resulting in degeneracies of 
F-point planewave eigenstates in the (d) nearly-free electron 
bandstructure along K — V — M axes. 


TABLE I. Character table of the F-point D^h group. Fi-g in 
parenthesis are corresponding IRs in the convention of Bethe 
notation, in which the extra IRs after the double group ex¬ 
tension are Fr-g not listed here. The basis functions, s and 
Px,y,z orbital configurations, and important invariants are also 
listed. The assignment of plus and minus superscripts to rep¬ 
resentations follows the convention of even and odd parity 
with respect to the operation of in-plane mirror reflection 
Ofi of C4. Superscripts on orbitals indicate relative sign be¬ 
tween orbitals on different (upper/lower) sublayers. Note that 
k = (kx, ky) = k{cos (f>, sin 0) 



Cl C2 C3 C4 C5 Ce 

basis 

orbitals 

invariants 

rr(ri) 

111111 

1 


k^ 

r2+(r2) 

11-111-1 

xy 


az, sinSc;/) 

r 3 +(r 6 ) 

2-10 2-10 

{x,y} 


{kx,ky} 

rr(r 4 ) 

1 1 1 -1 -1 -1 

xyz 



r 2 “(r 3 ) 

1 1-1-1-11 

z 

S~, pt 


r 3 “(r 5 ) 

2-10-21 0 

{xz,yz} 

{pZ,Py} 

{'XxjO'y} 


Here, similar to the IR labeling convention discussed 
above, even (odd) parity of orbital configuration within 
the unit cell upon application of au is denoted by the 
plus (minus) superscript. For example, p~ of F]^ consists 
of Pz orbitals of atoms in the two sublayers with opposite 
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wavefunction phase orientation. Because the atoms are 
from group-III and VI, these s and p orbitals dominate 
the low energy electronic structure covering the critical 
band gap region. Their limited combinations in Table |T] 
indicate that only four IRs (T^, T^, T;^ and T^) are 
relevant through out this paper. 

The right-most column of Table|T]lists several low order 
invariants according to their behavior under the twelve 
symmetry operations. For example, the function sin3(^ 
belongs to TJ, where (p is the polar angle of the 2D vector 
k with respect to the T—M axis. It transforms to — sin 3(j) 
under the six operations in classes C 3 and Cg that involve 
rotations along in-plane axes. These invariants are low- 
order terms in the expansion of the Hamiltonian based 
on k • p theory and will be used for the discussion of band 
dispersion near the zone center. 


B. Nearly free electron model 

In order to understand the origins of energy dispersion 
in the true band structure, we begin with a symmetry 
analysis of the nearly free electron (NFE) band structure, 
which we show along with corresponding symmetries of 
the zone center states in Fig. m- 

The lowest-energy state at the F-point (green hexa¬ 
gram) is of F)*" symmetry since its wavevector is the Bril- 
louin zone origin, labeled [00]. The first excited states 
(blue square) are sixfold degenerate, with their wavevec- 
tors corresponding to the six reciprocal lattice points 
nearest to the origin [Fig. [^c)] and denoted by [01]. 
By examining their characters under all the D^h symme¬ 
try operators, it is easily found that these first excited 
states are composed of 2 F 3 0 2F(‘. In fact, all reciprocal 
lattice points sitting on one of the three axes of in-plane 
180° rotation operators possess this symmetry. All other 
reciprocal lattice vectors share the same symmetry with 
the second excited states, denoted by red circles in Fig. 
m and labeled [ 11 ] (the wavevector sum of any two 
neighboring [01] states). These are also six-fold degener¬ 
ate and have the symmetry of 2r^ 0 F J 0 F]*". Given the 
fact that we are exclusively considering the low-energy 
bands, only the [ 00 ] and [ 01 ] states are of interest (the 
energy of [11] planewaves is nearly 30 eV). 

In the 2D NFE model, planewave eigenstates are al¬ 
ways even under the in-plane mirror reflection operator 
ah- In a physically real lattice where the configurations 
of atomic orbitals have even or odd parities with respect 
to ah, the Fg states of [ 01 ] can transform into Fg , while 
the F]'' states of both [ 01 ] and [ 00 ] can transform into F^ 
(see the sign of characters for € 3-5 that take z —>■ —z in 
Table |!|. 

The true crystal lattice potential modifies this 
planewave dispersion while maintaining the symmetry. 
First of all, the finite thickness of the real monolayer 
(analogous to confinement in a quantum well) results 
in a series of subbands that originate from the primor¬ 
dial 2D states. Secondly, the in-plane potential breaks 


the F-point degeneracy of different IRs in the NFE band 
structure, and mixes the eigenstates into a linear combi¬ 
nation of planewaves belonging to the same IR. Numer¬ 
ical procedures utilizing planewaves as basis functions, 
such as the empirical pseudopotential methocP^ and 
density-functional theory (DET) packages like Quantum 
ESPRESSO,l^ can be used to examine these planewave 
origins. 

Our last remark about this NFE model is that the 
momentum matrix elements between reciprocal lattice 
planewaves - even if they are allowed by symmetry 
- are nonvanishing only for degenerate symmetrized 
planewaves. For example. Table |l| indicates that the mo¬ 
mentum operator polar vector components and ity 
(we use this notation to avoid confusion with the Px,y,z 
orbitals) belong to the IR Fg , and their coupling between 
states belonging to F]*' and F]]" is apparently allowed by 
symmetry. Momentum matrix elements between states of 
these representations, if both are composed of [ 01 ] sym¬ 
metrized planewaves, clearly have magnitude on the or¬ 
der of 2'iTh/a (where a is the lattice constant). However, 
the coupling between F]^ of [ 00 ] and F]]" of [ 01 ] by nx,y 
vanishes due to the structure of their oscillatory wave- 
functions. We will show in the following subsection that 
this fact is essential in determining the band dispersion, 
because if two bands are coupled by the k-p perturbation, 
they tend to energetically repel each other. The cou¬ 
pling strength, and hence quantities like effective masses, 
thus sensitively depends on whether the vestigial origin 
of these bands is the same planewave state in the NFE 
model. 


C. Origin of valence band distortion, 
alias the ‘caldera’ 

The empirical tight binding band structure of GaSe in 
the vicinity of the band gap is shown in Fig. Besides 
the lowest F^ conduction band, five of the nine valence 
bands are plotted, including the non-degenerate highest 
valence band F]*" and the two pairs of doubly-degenerate 
Fg . It is interesting to note that the remaining four 
lowest valence bands not shown in Fig. [^are composed 
of two pairs of {r]'',r^} (see, for example. Fig. 1 of 
Ref. [H]). These states, differing by their definite parity 
with respect to ah, are each split by the sublayer bond¬ 
ing/antibonding energy. This heuristic can be confirmed 
by correlating the band splitting with the amplitude of 
P 2 orbital components (responsible for the strong inter¬ 
sublayer a bonding) in each pair. Likewise, the splitting 
of F^ states is determined by the tt bonding/antibonding 
energy difference between the even and odd configura¬ 
tions of purely Px and Py wavefunctions. 

Our first-principles calculation using norm-conserving 
pseudopotentials shows that the F]*" highest valence band, 
together with the remaining four nondegenerate lowest 
valence bands just mentioned, originate from the NFE 
[00] planewave. Their common origin can be understood 


4 


in terms of subband formation induced by electron con¬ 
finement to the quasi-2D atomic lattice with finite thick¬ 
ness, as discussed above in Sec. |IIB[ For example, in the 
F]^ highest valence state, this [ 00 ] component accounts 
for « 2/3 of the wavefunction amplitude, while only 30% 
comes from the [01] planewave. On the other hand, the 
two pairs of F^ valence bands, as well as the lowest 
conduction band F^, originate from the [ 01 ] planewave 
(more than 90% in Fg and 80% in F^). 

In light of these wavefunction compositions, we can 
conclude from the argument at the end of the previous 
subsection that symmetry-allowed momentum matrix el¬ 
ements coupling these bands, represented in Fig. [^by 


-Pi “ (^2 jflx.ylrg ), 

mo 

(red arrow) 

(1) 

P2 = —(r+l7r,,,ir+), 
mo 

(short pink arrow) 

(2) 


are quite different in magnitude. Specifically, the proba¬ 
bility associated with the former ([Pi J^) is about an order 
of magnitude larger than that of the latter (IP 2 P, ~ 0.3^ 
smaller). Therefore, the effective mass m* of F^ is much 
smaller than the free electron mass mg (similar to the sit¬ 
uation of the F-point conduction band minimum in many 
cubic systems such as GaAs and Ge), while m* of the F]'' 
valence band is slightly smaller than mg due to the weak 
upward repulsion from F;]" - but still positive! 

This unusual dispersion (which would otherwise result 
in a smaller, or closed, bandgap) is counteracted by the 
interplay of several secondary factors. First of all, down¬ 
ward repulsion comes from upper conduction bands with 
F;/ symmetry indicated by P 3 in Fig. [^[with similar def¬ 
inition as Eq. ([^]. Our DFT calculation shows that the 
third-lowest conduction bands, with F]/ symmetry, origi¬ 
nate from the [01] planewave. Close to the F-point, their 
downward repulsive influence on the F/" highest valence 
band (again, dominated by the [ 00 ] planewave) cannot 
compete with the upward contribution from the close- 
by F;/ valence band. Further away from the F-point, 
however, this upward repulsion quickly vanishes, because 
the F;]' lower valence bands are strongly repelled down¬ 
ward by the dominant [ 01 ] planewave component of the 
F[/ second-lowest conduction band (indicated by P 4 in 
Fig.[^. At large enough k, the F;/ conduction band wins 
over the valence band with the same symmetry, and over¬ 
whelms the (positive) quadratic free electron dispersion 
of the upper valence band. A large energy gap is thus 
opened. 

In summary, the ‘caldera’ shape of the upper valence 
band in the vicinity of the F-point is a result of the com¬ 
petition between the lower valence band and upper con¬ 
duction bands with the same F]]" symmetry, similar to the 
origin of the possible indirect bandgap in phosphorene 
along the zigzag direction.!^ The energy dispersion of the 
F[/ upper valence band can be analytically expressed via 
this k • p analysis and second-order perturbation theory 



FIG. 2. Spin-independent tight-binding bandstructure 
around the zone center along the K — V — M directions, in 
the vicinity of the bandgap. Here the j/-axis (energy) is con¬ 
tracted within the gap region for better illustration. IRs of 
zone center eigenstates are shown. The double-headed arrows 
indicate dominant k-p interactions responsible for a ‘caldera’- 
shaped highest valence band depicted by the 3D illustration 
inset. 

by 

^ ^ lP 2 pfc^ _ . 

2 mo Ea Eb ' ^ ’ 

where, in the second term, 

p _ T^VB _ (pVB _ \Pi?k'^ A 

is the energy difference between the F]/ highest valence 
band and the F;/ valence band. The dispersion of the 
latter is taken into account by the last term of Eq. Q, 
where P 4 has a definition similar to P 2 in Eq. ([^ but 
with an amplitude that rivals Pi. 

Note that in the large-fc limit, the second term of 
Eq. (§ approaches a constant, and its contribution to 
energy dispersion vanishes. In the last term of Eq. (|^, 
P 3 has similar amplitude to P 2 , and the energy denomi¬ 
nator 

Eb = {E^.f - Pk^) - EY^ (5) 

is the difference between the F[]" conduction bands and 
the r/" highest valence band. The —fik"^ < 0 term cap¬ 
tures the negative effective mass of the F;/ conduction 
bands due to k • p suppression from even further bands, 
and results in higher order dispersion of the F]/ valence 
band beyond Such higher order contributions, espe¬ 
cially the quartic term, are important for valence band 
edge hole states with relatively large wavevectors, e.g. at 
the caldera ‘rim’. 
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Given that Ea is almost an order of magnitude smaller 
than Eb at the F-point, it is very unlikely that the last 
term in Eq. ([^ would overcome the sum of the first and 
the second terms at the zone center and render an ordi¬ 
nary parabolic hole band with negative m*, as verified 
by both the empirical tight binding method and first- 
principles calculation. However, considering the possible 
deficiencies of numerical procedures (for example, under¬ 
estimation of the bandgap in DFT), we cannot unequiv¬ 
ocally assert the true nature of the valence band edge 
dispersion. Furthermore, as we will show in the follow¬ 
ing section, such dispersion reversal of the highest va¬ 
lence band at the zone center is sensitive to the spin- 
orbit coupling strength. Ultimately, the existence of this 
valence band ‘caldera’ must be empirically verified 
by experiment, such as optical spectroscopy (see Sec. 
to detect the divergent DOS associated with high-order 
dispersion, or angle-resolved photoemission spectroscopy 
(ARPES) to directly probe the shape of the valence band. 

With our fundamental understanding of the caldera’s 
origin, we can also easily explain its gradual disappear¬ 
ance and transition into an ordinary parabolic valence 
band edge when m onola yers are stacked and the system 
evolves toward bullPi^^. As the layer thickness increases, 
every band of the monolayer bandstructure develops a se¬ 
ries of subbands corresponding to different van der Waals 
bonding with all possible interlayer phase configurations. 
The splitting of these subbands is therefore independent 
of all k • p interactions, and rather strongly relies on the 
wavefunction amplitude of vertically-distant chalcogen 
anion orbital component. The T]^ valence band edge 
has a large anion orbital component!^ that diminishes 
away from the F-point,!^ tending to raise the bottom of 
the caldera more than the rim. Furthermore, and just 
as important, this subband splitting changes the energy 
denominators in Eq. causing weaker (stronger) re¬ 
pulsion from lower (upper) Fj^ bands. These altered in¬ 
teractions result in a gradually shallower caldera, which 
eventually disappears as bulk conditions are approached. 

Lastly, we comment on the weak anisotropy of the va¬ 
lence band, depicted in the inset of Fig. By tak¬ 
ing into account even higher-order contributions beyond 
the circular extrema given by Eq. (|^ at nonzero k 
can be modulated by a term proportional to fc® cos Gfj) 
(Ref. [H]). The resulting six valence band maxima occur 
on the T — K axes, whereas saddle points exist (nearly 
equidistant from the zone center) on the T — M axes. 


III. SYMMETRY AVEC SPIN 

A. Zeroth-order split-off states and spin mixing 

The inclusion of spin-orbit coupling into the Hamilto¬ 
nian can be treated within the framework of perturba¬ 
tion theory. Fig. depicts the important changes to 
the GaSe band structure among the five highest valence 
bands before (left) and after (right) spin-orbit interaction 


Dresselhaus 



FIG. 3. Evolution of valence bands before (left) and after 
(right) introducing spin-orbit interaction. Gray arrows indi¬ 
cate spin-orbit induced broken degeneracy of the F^ valence 
band. The symmetries of zone-center spin-dependent eigen¬ 
states are described by mixed basis functions given in Tablejl] 
Here, X = yz and Y = xz. 


is considered. The spin-independent bands are identical 
to those in Fig. except we have shortened the ener¬ 
getic distance between F)*" and F;)" for clearer illustration. 
Our incorporation of spin-orbit coupling into the tight- 
binding formalism follows the approach of Ref. m by 
considering on-site spin-orbit parameters of Se and Ga 
atoms. 

The most significant hallmark of spin-orbit coupling is 
the splitting of the F^ bands. Both of these bands have a 
double orbital degeneracy before including spin, similar 
to the top of the valence band in many cubic semicon¬ 
ductors in which the threefold degenerate valence band 
(sixfold if spin is included) is broken into the four-fold 
degenerate heavy- and light-hole bands and the doubly 
degenerate split-off hole bands. Here, taking F^ as an ex¬ 
ample, the two orbital states can be denoted by the axial 
vector component basis functions as {X = yz,Y = xz}, 
see Table They are coupled by the zeroth-order spin- 
orbit perturbation x p • ct, so that the split-off 

energy is determined by the off-diagonal matrix element 


A(F3-) 


2mgc2 


(A| 


/au d 

dV d\ 

\ dx dy 

dy dx ) 


in, 


( 6 ) 


and the spin-dependent eigenstates are a mixture of X 
and Y as shown in Fig. 3. The discussion for F(^ is the 
same, by replacing the basis functions {X, Y} in Eq. ^ 
with polar vector components {x,y}. Note that, due 
to the proximity of F;)" and Fg , this lowest-order spin- 
orbit interaction pushes down the lower pair of F^ and 
raises up the upper pair of Fj^ so much that their energies 
even switch order (crossing of the two gray arrows in the 
middle of Fig. |^. 

The lowest order spin-orbit interaction also couples 
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the F;^ to the F]^ highest valence band, since the 
k-independent spin-orbit invariants {ax,cry} belong to 
Fg = Fg 0 F]^. This much weaker non-degenerate per¬ 
turbation results in negligible energy shift compared with 
the split-off energy of Fg , but it does significantly alter 
the spin-dependent eigenstates by inducing spin-mixing. 
If we denote the spin-independent eigenstate of F^ by the 
scalar basis function 1 as shown in Fig. the unitless 
spin mixing coefficient a is 




a = 


4m^c^Ei+_3- 


-( 1 | 


rdv d 

dV d\ 

\ dy dz 

dz dy J 


1^)- (7) 


Here, £'i+_ 3 - is the energy difference between F^ and 
Fg , and their proximity (only several hundred meV) 
leads to a relatively large value of a. Our tight-binding 
calculation shows that for F^, the spin-mixing defined by 
the total square amplitude of the minority-spin compo¬ 
nents is « 8 % at the F-point, falling only slightly to 
~ 5% at the caldera rim. For comparison, we also have an 
interaction between the F|^ valence bands and F^ lowest 
conduction band that induces spin-mixing, yet it is much 
less pronounced due to the much larger energy denomi¬ 
nator; our calculation gives 0.1% . Sin ce the strength of 
Elliott-Yafet (EY) spin relaxatiorP^EI] jg governed by the 
spin-mixing amplitude, the resulting EY relaxation rate 
for spin-polarized conduction electrons is nearly two or¬ 
ders of magnitude smaller than holes near the top of the 
valence band. Another important spin relaxation mech¬ 
anism, Dyakonov-Perel (DP), is due to spin splitting of 
every band and is discussed in the next subsection (Sec. 

hIb ). 

Due to the mixing of X and Y components into the 
F)^ highest valence band by spin-orbit coupling, opti¬ 
cal selection rules allow in-plane polarized electromag¬ 
netic radiation to connect states across the fundamental 
bandgap. This mixing is therefore especially important 
for experiments, since it allows carrier generation with 
normally-incident band-edge illumination.EIlll We elab¬ 
orate on radiative transitions and corresponding selection 
rules relevant for conduction-band spin polarization via 
optical orientation in Sec. |V] 

Lastly, we note that in Fig. the zeroth-order spin- 
orbit coupling reduces the depth of the ‘caldera’ distor¬ 
tion of F^*" highest valence band. As discussed in the pre¬ 
vious section, this distortion is due to repulsive competi¬ 
tion between the valence and conduction F|^ bands. Since 
the energy differences to F)*" are modified by the split- 
off energy A(F(^), dispersion close to the Brillouin zone 
center is sensitive to spin-orbit strength. Moreover, the 
spin mixing of F;^ components into the highest valence 
band allows its coupling to bands with F^ symmetry via 
the k • p perturbation, incurring an additional downward 
repulsion contribution from the conduction band mini¬ 
mum. In fact, for the tight-binding GaSe model we use, 
positive effective mass at the F-point disappears for only 
«35% stronger spin-orbit strength. With this concern, 
we reiterate our lack of absolute certainty on the ques¬ 
tion of whether the quantitative nature of the highest 


valence band is a ‘caldera’ shape or an ordinary hole-like 
paraboloid. 


B. fc® spin splitting: Dresselhaus effect 

A remaining feature of the spin-dependent band struc¬ 
ture shown in Fig. is the splitting of opposite spin sub¬ 
bands for k along F — K{K') directions. Although this 
splitting vanishes on the T — M axes (the C 2 v group of 
the wavevector for points on these axes only has two- 
dimensional double group IRs), it generally persists for 
an arbitrary choice of k. This splitting is similar to Dres¬ 
selhaus spin-splitting in zincblende semiconductors,^ 
where it vanishes along F — Likewise, it has the 

same fundamental cause, namely the absence of space 
inversion symmetry that allows a spin-orbit-induced ef¬ 
fective internal magnetic field in reciprocal space. In this 
subsection, we analyze the symmetry of this higher-order 
spin-orbit effect using the method of invariants.!^ 

There are two potential approaches within this frame¬ 
work. One could extend single group to double group 
notation, in which the Pauli matrices of spin orbit in¬ 
teraction are embedded as invariant matrices, while k- 
dependent terms are filled into the correct position in 
the matrices as invariant components.!^ However, for 
simplicity we follow an alternative approach, considering 
both k and Pauli matrices (and their combinations) as 
invariant components within the single group notation.!^ 

We begin by taking the F)^ highest valence band as an 
example. The direct product of this IR with itself gives 
F)^ GF^ = F)*", indicating that any term including ax or 
Oy is forbidden. These two Pauli matrices behave as axial 
vectors (see Table |l]), which are odd under the operation 
of in-plane mirror reflection au- The same is true for 
their combination with kx and ky (both are even under 
ah)- This existence of the in-plane mirror reflection is 
a critical factor that determines the form of Dresselhaus 
spin splitting. In wurtzite systems sharing the threefold 
rotation symmetry but lacking an as a group element, 
ax or ay are allowed in the spin-dependent Hamiltonian 
and the lowest order Dresselhaus field is linear in kx and 



The remaining Pauli matrix, tJz, belongs to FJ. Its 
combination with polynomials of kx = k cos (j) and ky = 
ksin(j) that also belong to F^ are symmetry-allowed in¬ 
variant components for the F)*" band (F^ G F^ = F)^). 
We find that sin 3(j) = 3 cos^ (j) sin (j) — sin^ (j) (requiring 
k^ dependence) is the lowest-order invariant component 
that belongs to F 2 . With all these concerns, we conclude 
that the lowest-order Dresselhaus term can be written as 

77so = 7 i sin 30 , ( 8 ) 

where 71 is the spin-orbit coupling strength coefficient 
for this band. The form of Eq. Q is universal for 
all other bands (with different 7 coefficients), including 
those originating from the two-dimensional Fg bands. 
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due to the constraint imposed by ah and three fold rota¬ 
tional symmetries. Fig.|^a) is a schematic representation 
of Eq. , where blue arrows represent the effective inter¬ 
nal magnetic field. The field magnitude is proportional 
to and alternates direction between parallel and anti¬ 
parallel relative to the surface normal as a function of 
polar angle. 

The spin-splitting parameter 71 in Eq. (|^ can be cal¬ 
culated using fourth-order perturbation theory. As in the 
corresponding calculation for III-V semiconductors, this 
quantity involves three matrix elements of the k • p inter¬ 
action and one of the spin-orbit coupling 47^72 VE x p-tr. 
Alternatively, spin-orbit can be taken into account to all 
orders by using the exact split-off energies in third-order 
perturbation term denominators.^ 

As discussed in the previous section, the k • p terms be¬ 
long to Tg and can couple the T^^ highest valence band 
only to bands with Ej^ symmetry. The Tj^ intermediate 
states, in turn, can also be coupled to each other by k-p, 
since Tj^ 0 E;^ = E^ 0 E^ 0 E;^. The dominant third- 
order k ■ p perturbation paths are therefore those shown 
in Fig. mb), and include the Ej^^ valence band and the 
lowest rjc conduction band (the same as those given in 
Fig.l) as intermediate states. For simplicity, their ener¬ 
gies with respect to E)*), in ascending order are denoted W 
Ai_ 4; clearly, E2 — Ei is equal to A(E3 ) given in Fig. ^ 
and A4 — A3 is the split-off energy of the conduction E^ 
states (with a similar definition). By summing over all 
paths in Fig.j^b), one obtains 


7i 


ml 



{l\TTy\i){i\TTy\j){j\TTy\l) 

EiEj 


\P2Qp^\ 


\E2E3 



( 9 ) 


where P2 and P3 are given in Fig. and Eq. ( 1 ^, and 
Q is the k • p matrix element Due 

to the partitioning of x and y components in spin-orbit- 
split Ej^ wavefunctions (as shown in Fig. Q couples 
energy levels only between Ei and A4, and between E2 
and A3,1 ^so that in Eq. (|^ there are no terms with en¬ 
ergy denominator E1E3 nor E2E4. Note that the corre¬ 
sponding calculation for the E^ lowest conduction band 
involves dominant contribution from the Eg^ ^ bands, and 
is similarly dependent on the spin-orbit-induced splitting 
A(rg-). 

The out-of-plane orientation of the internal magnetic 
field in group-III metal-monochalcogenides resembles the 
case of IIEV zincblende semiconductor [ 110 ] quantum 
wells, except that in the latter system Dresselhaus split¬ 
ting is linear in k (see Ref. [ 60 ]). In that case, spins 
oriented normal to the plane have relatively much longer 
lifetime than in-plane spins because they are eigenstates 
in the Dresselhaus field, and thus the Dyakonov-Perel 
(DP) spin relaxation mechanism vanishes.!^ We there¬ 
fore expect strong spin relaxation anisotropy in the 
monochalcogenides for the same reason. Specifically, 



FIG. 4. (a) Spin-orbit-induced fc-cubic Dresselhaus splitting 
[Eq. of the spin-np and spin-down subbands is represented 
by the red manifold. The green hexagonal plane is similar to 
the reduced Brillouin zone with F — A and F — M axes shown. 
Blue vectors represent the Dresselhaus internal magnetic field. 
Note that the out-of-plane field direction is reversed for K and 
K' and vanishes along F — M. (b) Dominant perturbation 
paths (the two dashed loo^ for the calculation of the coef¬ 
ficient 71 in Eqs. (|^ and Q. Horizontal lines represent the 
spin-dependent F-point states, with their single group origin 
listed on the right. E 1-4 are the energy values of the spin-split 
states with Fg" single group origin, with respect to the energy 
of F^'j,, in ascending order. Pg and P 3 are the same as they are 
in Fig. [^ and Q is the momentum matrix element between 
conduction and valence bands with the same Fg" symmetry. 


when the spin orientation is chosen to be in-plane, spin 
relaxation is governed by DP and proportional to 71. On 
the other hand, spin with out-of-plane orientation is not 
subject to precession in the Dresselhaus field and DP is 
absent; in that case, EY (determined by the spin-mixing 
coefficient a) then dominates spin relaxation. 

Another feature of the [ 110 ] zincblende quantum well, 
enabled by the k-linear spin splitting and fc-quadratic 
band dispersion, is the possibility to generate a so-called 
‘persistent spin helix’. E21E1 This static spin texture re¬ 
sults from the cancellation of ^-dependence in the ratio 
of spin-orbit field to group velocity determining preces¬ 
sion angle. The monochalcogenide hole states at the E]'' 
valence band maximum have a large dispersion com¬ 
ponent, and together with Dresselhaus splitting thus 
can also support a ‘spin helix’. However, the underlying 
threefold rotation symmetry eliminates its ‘persistence’ 
in the presence of scattering. 

As long as the threefold rotational symmetry is cap¬ 
tured, terms higher than cubic in k are allowed, if con¬ 
structed by combining Eq. (|^ with an additional invari¬ 
ant component belonging to Ej''. For example, noting 
that belongs to Ej*", the possible presence of fifth-order, 
seventh-order, etc. in k cannot be ruled out. These 
higher-order contributions become prominent at large fc, 
especially for bandedge states on the rim of the valence 
band ‘caldera’. If such higher order terms have the op¬ 
posite sign as the lowest-order cubic term, the Dressel¬ 
haus spin splitting would appear sub-cubic as k increases. 
However, we must emphasize that in this system, a linear 
Dresselhaus term is forbidden by symmetry. Any con¬ 
trived attempt to fit a numerically-calculated spin split- 












ting using only linear and cubic terms, as from DFT in 
Ref. [ 65 ], is physically unjustified and susceptible to mis¬ 
leading evidence. Mistaken conclusions are an inevitable 
consequence of artificial symmetry breaking stemming 
from e.g. numerical rounding errors (associated with un¬ 
avoidably representing the irrational values of atomic po¬ 
sition components of the hexagonal Bravais lattice in a 
finite-digit scheme) at some arbitrary level of precision. 


IV. ORBITAL MAGNETISM AND EFFECTIVE 
g-FACTOR 

As previously mentioned, the spin-orbit-induced split¬ 
ting of lower Fg valence bands is reminiscent of the 
r4 —^ r70r8 valence band splitting in bulk cubic semi¬ 
conductors. In that familiar case, the valence-band split¬ 
ting allows a nonzero orbital contribution to the mag¬ 
netic moment and a conduction-band g-factor substan¬ 
tially different from the spin-only value of 2 , with smaller 
band gap and larger split-off energy strongly enhancing 
the correction F 1 1 66 | 67 ] jjeuristically, we can surmise that 
large orbital g-factor in a given band follows from cou¬ 
pling to nearby bands which originate from states whose 
initial orbital degeneracy has been split by spin-orbit in¬ 
teraction. In MX, the lower F^ valence bands are indeed 
split by an amount on the order of the band separation 
to the upper Fj*" valence band, suggesting a large orbital 
g-factor there. To test this assumption, we therefore un¬ 
dertake the appropriate calculation below. 

In a magnetic field, we make the Peierls substitution 
by writing the canonical momentum vector tt —tt -|- eA. 
An out-of-plane field i?^z has vector potential A = B^xy 
(in the Landau gauge) so that the Hamiltonian term giv¬ 
ing rise to orbital diamagnetism is The 

expectation value of the operator in parenthesis can be 
evaluated in the Bloch band basis to yield the orbital 
magnetic moment (in units of |^ = ^) that competes 
with the spin magnetic moment to give an overall effec¬ 
tive g-factor < +2. 

Unfortunately, the position operator x cannot be eval¬ 
uated directly due to the indeterminate nature of matrix- 
element integrals over delocalized and non-normalizable 
Bloch waves that extend to infinity.l^ Instead, using the 
Ehrenfest theorem (tTj,) = mvx = = ^[XjH] and 

taking matrix elements in a band basis \'4’n), we obtain 
(V'n|i|^/’n') = in R-of. jSS] . Thus, the or¬ 

bital “g-factors” of each band can be evaluated using the 
operator with matrix elements 


{n\gorbit\n'') 


4i 

m ^ En- En' 

n'^n 


In a symmetric gauge, iTyX —>■ ^{TTyX — Tr^y), with appro¬ 
priate change to the numerator in Eq. We note that 
in the tight-binding formalism it is straightforward to 
use momentum operators constructed directly from the 



FIG. 5. K — V — M orbital magnetic moment of both spin 
states in the Fj*" valence band. Units chosen indicate orbital 
contribution to effective g-factor. Inset: Orbital magnetic 
moment of uppermost valence band within the full Brillouin 
zone, highlighting KjK' valley magnetism. 


Hamiltonian "H, via tt = although this ignores 

intra-atomic orbital overlap contributions while main¬ 
taining gauge invarianc^^^HZH 

Immediately we see that the coupling to bands gov¬ 
erning the orbital g-factor is through matrix elements of 
orthogonal components of momentum. [An interesting 
example is phosphorene,^ where due to strongly broken 
symmetry all orbital bands are nondegenerate, and fur¬ 
thermore along the high-symmetry axes in fc-space the 
momentum matrix elements of tTx and Hy are never si¬ 
multaneously nonzero, so that orbital g-factor of every 
band vanishes.] According to Table |lj the polar vector 
components transform like F]]", so coupling to the 
Fj^ upper valence band is dominated by the Fj]" lower 
valence band (and F^ upper conduction bands, but to a 
lesser degree due to the large energy denominator). 

As discussed in Sec. mE the spin-orbit field vanishes 
along T — M and spin degeneracy is preserved. For this 
axis, orbital magnetic moment for spin up and spin down 
in the Fj'' upper valence band, calculated using first-order 
degenerate perturbation theory to evaluate Eq. ( 10 ), are 
equal in magnitude but opposite in sign as expected (see 
Fig. §. However, lattice inversion a symrn etry allows 
valley-dependent magnetic interactions, F^F3 | jg 

manifest in the behavior of orbital magnetic moment 
along the T — K axis, where spin degeneracy is lifted. 
Instead of sharing equal magnitude orbital magnetic mo¬ 
ment with a nearly degenerate state at the same k, here 
the opposite magnetic moment is found with its time- 
reversed partner (along T — K'). Thus, all three K points 
(oriented at 120 ° in fc-space) have identical magnetic mo¬ 
ment (and at jgj « 3 , several times larger than at F), with 
all three complementary K' points having the opposite 
moment, nearly independent of spin (see inset to Fig. 
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for g-factor in the full BZ). The lowest conduction band 
shows a similar dichotomy between behavior of magnetic 
moment along T — M and V — K, but the orbital g-factor 
for electrons is quite small, reaching only approximately 
\{gorbit) \ =0.3 (not shown). 

This relatively small orbital correction to effective g- 
factor is contrary to our initial heuristic expectation of 
a large correction in the T)*" band based on spin-orbit 
splitting of the nearby T;^ band. However, while the de¬ 
nominators in Eq. are indeed small and unequal as 
expected, the momentum matrix elements between these 
two bands is weak. This can be understood by recalling 
from Sec. IIB that the dominant planewave origin of the 
T)^ upper valence band is from [00], whereas the lower T^ 
valence bands are rooted instead in the [01] planewaves 
[see Fig. [^d)]. Momentum matrix elements between 
these two bands, dominated by contributions from dif¬ 
ferent planewave origins, are inherently weak. 


V. OPTICAL ORIENTATION 
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FIG. 6. Spin-resolved optical absorption of right-handed 
circnlar polarized light resniting in excitation of electrons 
into the conduction band, and resulting spin polariza¬ 
tion (black line, right axis). Critical band edge energies are 
indicated by dashed or dotted lines. 


With appropriate selection rules and spin-orbit split¬ 
ting of the otherwise degenerate valence band, conduc¬ 
tion electron spin polarization can be generated via in¬ 
terband dipole excitation.l^ This “optical orientation” 
process is most clear in cubic III-V semiconductors with 
a (£ = 1) p-like valence band, where one j = 1/2 band 
is split off from the remaining j = 3/2 states by low¬ 
ering an energy A. Then, circularly-polarized photons 
with energy hui at the bandgap Eg can optically excite 
electrons into the conduction band, with matrix element 
asymmetry for spin up to spin down at the zone center 
of 3:1. This ratio is determined by the orbital wave- 
function components (e.g. Clebsch-Gordan coefficients), 
giving spin polarization P = = 50% for fno = Eg¥^ 

Higher spin polarization can be obtained in 2D epitax¬ 
ial quantum wells made from the same materials, where 
confineme nt ind uces splitting of the light and and heavy 
hole states.l^^^^ 

Optical orientation in bulk GaSe was explored both 
theoreticalljli^ and experimentalljiA^ in the 1970s (and 
more recently in Ref. [3H])? but the symmetry of single¬ 
layer GaSe is expected to yield different results. 

A. Absorption Spectroscopy 

We calculate the corresponding optical orientation in 
the conduction band of monolayer GaSe within the tight- 
binding framework, by sampling wavevectors within the 
full Brillouin zone and compiling the squared (circularly- 
polarized) optical matrix elements in 

A^/^{huj) OC - {Ec - En)), (11) 

n^c 

where 'k± = tt^, ± iffy determines the chirality of the 
circularly-polarized electromagnetic field; this approach 


restricts our analysis to direct transitions only. We sub¬ 
sequently calculate the conduction spin polarization via 

p _ A^-A-^ 

— Af+Ai ■ 

The results from a sample of 10® pseudo-random k 
points is shown in Fig. where transition probability 
via right-handed circularly polarized light into the spin 
up (down) conduction band is shown in red (blue), to¬ 
gether with the spin polarization Pc in black. One strik¬ 
ing feature of these curves is the abrupt step-like absorp¬ 
tion edges due to the 2D DOS singularity caused by the 
‘caldera’ valence dispersion. Because of this behavior, we 
are justified in displaying the least-squares fit to the spec¬ 
troscopy with a piecewise quadratic function. Note that 
the divergent oc (—density of states associated 
with the quartic oc dispersion in the T/" valence band 
causes the near-band-edge absorption to initially fall as 
photon energy huj increases. The existence of the puta¬ 
tive ‘caldera’ valence dispersion can therefore be straight¬ 
forwardly supported by this one experimental signature, 
if it can be isolated from excitonic modification to the 
absorption spectrum. 

Due to dominant r/" character of the upper valence 
band, optical absorption of normally-incident circularly 
polarized photons at the fundamental band edge is sup¬ 
pressed (i.e. (Fj^lTg |r 2 ) = 0). However, the small quan¬ 
tity of -^{X E iY) spatial character with opposite spin 

mixed into the wavefunction from the T;/ lower valence 
band (see Fig. results in perfect -100% optical orien¬ 
tation, despite weak optical efficiency. 

Once the photon energy Klo is sufficient to excite elec¬ 
trons from the Tg band at approximately 3.35 eV, the 
conduction band spin polarization nearly fully reverses 
due to the much stronger optical matrix element from a 
dominant ;^(A1 ± lY) character. However, further in- 
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creases in energy by the split-off increment A(r;^) popu¬ 
late both spin states nearly equally, destroying the opti¬ 
cal orientation for hut >3.7 eV. Again, this transition is 
much more abrupt than in the case of cubic zincblende 
direct-gap semiconductors because of the higher zone- 
center density of states in each band. 


B. Valence band polarization via relaxation 


An intriguing phenomenon is enabled by the spin struc¬ 
ture of electron states and corresponding selection rules 
in these monochalcogenides. Refer to Fig. [Tja): Opti¬ 
cal excitation with normally-incident circularly polarized 
light of sufficient energy generates spin polarized carri¬ 
ers in the F^ conduction and the upper branch of the 
spin-split Fg lower valence band. Direct transition of 
the F;^ hole to the F]^ upper valence band is suppressed 
by the optical selection rules involving dominant char¬ 
acter of the wavefunctions of each band, so low-energy 
but large-momentum phonon emission is necessary for 
relaxation. This scattering event takes an initially spin- 
polarized hole near the F-point to states far away from 
the zone center, from where it can relax to the valence 
band maximum via a cascade of optical phonon emi ssion 
events, shedding tens of meV energy at each step.^^^*^ 
Due to strong spin-orbit-induced spin mixing (with am¬ 
plitude a, see Fig. in this band, the Elliott spin scat¬ 
tering mechanism results in virtually complete depolar¬ 
ization upon thermalization. However, the F^ conduc¬ 
tion band and F^*" valence band are connected by a strong 
TTj, optical matrix element with F^ symmetry, so radia¬ 
tive relaxation proceeds via linearly-polarized lumines¬ 
cence emitted in an in-plane direction. Since this direct 
interband transition is spin-conserving, oppositely polar¬ 
ized holes will remain in the F^*" valence band. 

This dynamic upper valence band hole polarization 
can be modeled via a rate-equation approach, similar 
to the simulation of recent time-resolve d lum inescence 
studies in thin but still 3D GaSe flakes.^^^^^ A simpli¬ 
fied equivalent three-level system is shown in Fig. [^a), 
where spin-conserving interband transitions (Fg o F^ 
optical orientation with generation rate constant and 
spontaneous relaxation rate constant R± , and F^*" —>■ F^ 
bandgap radiative relaxation with rate constant R^) ap¬ 
pear in solid arrows, while spin-mixing intraband spin 
relaxation (with timescales Tc^v,i for conduction, valence 
and lower-valence bands, respectively, indicated by sub¬ 
scripts) appears in dashed arrows. The phonon-assisted 
Fg —>■ F(^ inter-valence-band transition (with time con¬ 
stant T/i) appears as a solid horizontal arrow. 

We can write the following rate equations for the den- 

"T" /-i- 

sities of conduction-band electrons ric 


dnt 

dt 


dn^ 

dt 


=X^ 


nl - nj 

Tc 

nj - n\ 
Tc 


RzuIpI - R±nlpl, 
Rznipi - R±'nip\, 


( 12 ) 

(13) 




FIG. 7. (a) Schematic illustration showing optical orientation 
of spin-polarized conduction electrons via circularly-polarized 
electromagnetic excitation from the lower valence band and 
subsequent relaxation dynamics resulting in polarization of 
the upper valence band, (b) Spin-dependent carrier density 
in an equivalent three-level system modeled by Eqs. mE 
and (c) the corresponding spin polarizations. In normalized 
units, parameters are Tc = 80, r„ = 1, re = 1, Rz = 10, 
Th = 1, Dc = 5, Dt = 5, and = 0. 
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and (lower) valence band holes pV^- 

dpi 

dt 
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(14) 

(15) 

(16) 
(17) 


where the generation terms ^{Di—p[ ^){Dc— 

'n}J^) incorporate transition dependence on occupancy 
within a finite density of states Dc and D^. These non¬ 
linear equations are solved under steady-state conditions 
(all ^ = 0) with the added constraint of global neutrality 
dictated by conservation of charge, 

{nl + ni)-{pl+pi+pl+pi)=0. (18) 

Example results for steady-state carrier densities are 
shown in Fig. Bb) as functions of excitation inten¬ 
sity oc G^, where rounded relative parameter values are 
chosen for illustrative purposes, and perfect optical ori¬ 
entation (G'^ = 0) is assumed. Because direct (elec¬ 
tromagnetic radiative) transitions are faster than those 
that involve emission or absorption of large -momentum 
phonons during hole energy relaxation from to T)'', 
RzDc is chosen to be more than an order of magnitude 
larger than the inter-valence band momentum relaxation 
rate I/th- Furthermore, as discussed in Sec. 


Ill A 


our 

tight-binding calculations show that the F^ conduction 
band wavefunction has only 0.1% minority spin-mixing 
probability compared to the valence band at 8%, 
which is reflected in the relative values for the spin life¬ 
times Tc,v,l- 

It can be seen that at low excitation rates, the conduc¬ 
tion and upper valence band carrier concentration are 
proportional to VG7, whereas the lower valence band 
has a linear dependence. The former behavior is due to 
mutual Pauli exclusion suppressing recombination from 
these bands, whereas hole transitions from lower to up¬ 
per valence band are controlled only by the initial carrier 
density. 

With increasing excitation intensity, the carrier densi¬ 
ties saturate. In this regime, the conduction band spin 
polarization P^. approaches 100%, and the upper valence 
band spin polarization Py deviates substantially from 
zero, as can be seen in Fig. [^c). Because n| —>■ Dy 
as G^ —)■ oo when G'^ = 0, an exact analytic expression 
for the saturation polarization can be obtained by solu¬ 
tion to a polynomial equation of only quadratic order, 
but its utility to provide physical insight is limited by 
(grotesque) complexity. Nevertheless, under the assump¬ 
tion that Tc ^ Ty,h and DyR^Ty ^ I, we can simplify it 
to lowest order 


Using parameters from the simulation shown in Fig. 
Eq. (19) yields upper valence band polarization of -92%, 
in excellent agreement to the calculated « —92.45%. Un¬ 
der the same assumptions, the lower valence band polar¬ 
ization P^ « 2 Th+Tt again closely matching our 

numerical result of «-I-33.83%. 


VI. SUMMARY 

Our study of 2D MX monolayers has revealed the 
fundamental origins behind several properties of its elec¬ 
tronic structure. Throughout this group-theory analysis 
begun in Sec. m symmetry-allowed momentum ma¬ 
trix elements have played an essential role in determining 
e.g. the k • p coupling leading to valence band ‘caldera’, 
spin-orbit effect on band structure, orbital contribution 
to g-factor, optical absorption selection rules, and degree 
of carrier spin polarization from optical orientation with 
circularly-polarized electromagnetic radiation. We com¬ 
ment on the significance of each below. 

As we saw in Sec. m the valence band caldera is 
formed through a competition between repulsive interac¬ 
tions from bands at higher and lower energies. Impor¬ 
tantly, we found in Sec. |IIB| that since the eigenstates 
are confined to 2D, subbands originating from the pri¬ 
mordial NFE [00] planewave can produce the F]'' highest 
valence band, above the [01]-rooted Fg bands. The weak 
interband dipole interaction with these lower bands, the 
vanishing effect of the F^ conduction band just above, 
and interactions from even higher bands ultimately re¬ 
sults in the flat caldera. Two-dimensional confinement 
is the essential ingredient here, so this phenomenon is 
not expected to be unique to MX. In fact, although its 
symmetry is quite different from MX, phosphorene has 
an ultra-flat valence band along the zigzag direction for 
similar reasons.!^ 

Spin-orbit coupling induces several important features 
in the electronic structure. First of all, as detailed in Sec. 
|IIIB[ broken lattice inversion symmetry allows Dressel- 
haus spin splitting, and the remaining point-group sym¬ 
metry excludes any in-plane effective field components 
and splitting dependence linear in k. More importantly, 
as discussed in Sec. |III A[ the most obvious consequences 
of spin-orbit coupling are the broken twofold orbital de¬ 
generacy of the F^ bands, and spin mixing into the band- 
edge states. The latter effect determines the scattering- 
induced spin relaxation rate, and the former enables opti¬ 
cal orientation and an orbital contribution to the effective 
g-factor. 

The inequivalence of g-factor for states in the Brillouin 
zone toward K and K', discussed in Sec. IV has the same 
root as valley-spin coupling in TMDCs. There, however, 
the valence band maxima occur not at the F-point, but at 
the K and K' points per se. Although we do indeed see 
a weak valley-dependent g-factor in MA, the six shallow 
valence band maxima on the ‘caldera’ rim are invariably 
strongly coupled through intervalley scattering, making 
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experimental confirmation through spin Hall effect or val¬ 
ley Zeeman effect especially challenging. 

Because they determine the electromagnetic dipole 
symmetry, momentum matrix elements dictate the op¬ 
tical selection rules. We saw in Sec. [V] that in-plane 
momentum matrix elements connecting the and F^ 
band-edge states are inherently weak and only due to 
spin-orbit-induced band mixing of F;^ character into the 
upper valence band. Absorption of normally-incident 
transverse-polarized electromagnetic plane waves with 
huj at the bandgap energy is thus suppressed. This con¬ 
straint on optical absorption does not exist in the three- 
dimensional bulk counterpart material, so that thinning 
to monolayer drastically attenuates optical absorption. 
It is interesting to contrast this situation to TMDCs, 
whose transition from indirect-gap in bulk to direct-gap 
in monolayer results in a massive increase in optical 
absorption.l^ 

These optical selection rules can be exploited to gener¬ 
ate spin-polarized electrons in the conduction band with 


circularly-polarized photons, as discussed in Sec. V A 


When carriers are excited via transitions from the F^^ 
lower valence band, relaxation down to the upper valence 
band (by emission of electromagnetic radiation polarized 
normal to the plane and propagating in-plane for elec¬ 
trons, and via non-spin-conserving cascade mediated by 
phonon emission for holes) results in oppositely-pol arize d 


VB 


F)^ upper-valence band holes, as discussed in Sec. 

We note that our application of group theory to ana¬ 
lyze the symmetry properties of the MX monolayer can 
be straightforwardly extended to the case of stacked few- 
layer 2D flakes. Depending on the stacking configura¬ 
tion and number of layers, the inversion asymmetry of 
the monolayer (endowing many of the important spin- 
dependent qualities discussed in this paper) may be re¬ 
tained. Otherwise, time-reversal symmetry will preserve 


spin degeneracy and the Dresselhaus splitting will disap¬ 
pear. 

Other potential extensions of our theoretical frame¬ 
work include an incorporation of the mechanical lattice 
dynamics to determine phonon symmetries that could be 
used to expose the qualitative and quantitative aspects 
of various momentum and spin relaxation mechanisms 
relevant to our work. For example, the definite parity 
of all IRs with respect to reflection symmetry (Th decou¬ 
ples flexural phonons (odd parity) from participating in 
spin-conserving momentum scatterii^ whereas their ef¬ 
fect on EY spin-flip may be strong^S. The incorpora¬ 
tion of extrinsic invariants into our theory, such as mag¬ 
netic field, electric field, and strain, is straightforward 
once their transformation properties and corresponding 
IRs are identified. In particular, these external fields 
may break the in-plane mirror reflection symmetry of 
the monolayer lattice and induce linear Bychkov-Rashba 
spin-orbit couplin^^. The resulting interplay with the 
cubic Dresselhaus field may provide a potential scheme 
for electronically-actuated spin manipulation. 
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